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The elastic constants of a discotic lyotropic nematic liquid crystal are calculated by means of a 
pseudo-molecular approach as functions of the micellar shape anisotropy. By assuming that the 
temperature dependence of the ratio of the elastic constants comes from the temperature depen- 
dence of the micellar shape anisotropy, the theoretical predictions are connected with experimental 
measurements for the ratio K33/K11. This procedure permits to determine, in a phenomenologi- 
cal way, the temperature dependence for the ratio of elastic constants and for the micellar shape 
anisotropy near the nematic-isotropic transition in agreement with the experimental data. 

Lyotropic nematic liquid crystals (NLC) are formed by mixtures of amphiphilic moleculespj and a solvent (usually 
water), under convenient temperature and concentration conditions. The basic units of these systems are anisotropic 
micelles 2J. Under concentration change, or temperature change, the system can exhibit two uniaxial and one biaxial 
phases 3] . In the temperature-concentration phase diagram, the biaxial phase appears to be an intermediate phase 
between the uniaxial ones. The uniaxial nematic phases are known as discotic (Nj), and calamitic (N c ), depending on 
their magnetic anisotropy 2j. For what concerns the elastic properties of a NLC, one of the fundamental aspects of the 
research is the determination of the Frank elastic constants. For a distorted NLC the spatial variations of the director 
are small over molecular dimensions and a continuum elastic theory can be applied. In the bulk of a NLC, in general, 
can occur three basic types of distortions: splay, twist, and bend. To these distortions are associated, respectively, 
the elastic constants K\\, K22 and K33. These elastic constants have been the object of extensive experimental 
and theoretical studies specially when dealing with thermotropic NLC. An important parameter in this study, in 
compounds with a rigid molecular structure, is related to the molecular length-to- width ratio 0. 

Experimental results concerning lyotropic systemsQ show that the shape anisotropy of the micelles varies according 
to the composition and temperature. The micellar shape anisotropy (R) is defined as the ratio between the micellar 
diameter and the bilayer thickness. Its variation in a discotic nematic phase as a function of the temperature was 
determined in order to investigate a similar behavior of the orientational order parameter [(J. In this direction, it can 
be very instructive to investigate the dependence of the elastic constants on the anisotropy in the micellar shape of a 
lyotropic medium. In the literature Q, very few investigations have been reported on the elastic constants associated 
to the lyotropic systems. In particular, Haven, Armitage and SaupeQ have carried out precise measurements on 
the ratio of the elastic constants -K33 and K\\ as a function of the temperature for the discotic mixture of the 
decylammonium chloride (DAC'l), obtaining values comparable to those determined in nematic thermotropics. 

^From the theoretical point of view, it has been recently utilized a molecular approach to evaluate the micellar 
shape anisotropy dependence of the bulk elastic constants in a discotic nematic, by using the Maier-Saupe interaction 
law|9l Ho| . In that approach, the value for the ratio of the elastic constants K^/K\i was determined as an increasing 
function of the parameter R. The tendency and the order of magnitude of the ratio of elastic constants was found 
to be in a reasonable qualitative agreement with some available experimental results and with some calculations for 
thermotropics [Tlj . However, in that calculation the value for this ratio was limited to a saturation value K33/K11 = 2 
, thus preventing a broad comparison with other experimental data. 

In this paper we apply the molecular approach to calculate the elastic constants as a function of the micellar shape 
anisotropy having in mind a lyotropic nematic phase whose micelles are in the shape of a disk. We consider a sp ecial 
kind of interaction law which is a mixing between the Maier-Saupe and Nehring-Saupe interaction laws |l2t Il3| . For 
this interaction law the ratios between the elastic constants are not limited to these relatively small values for lyotropic 
NLC. Furthermore, we use the experimental data of ref. @ to obtain a phenomenological dependence of the ratio of 
the elastic constants as a function of the temperature, by assuming that this dependence comes from the temperature 
dependence of the micellar shape anisotropy. 

To determine the elastic constants of a discotic lyotropic we employ the pseudo-molecular approach, which is an 
approximate technique to determine the macroscopic properties of the system from the intermolecular interaction 
giving rise to the nematic phase^[0,^|. In this approximated method the elastic energy density can be written as 

/ = fa + Likrii.k + NijkTiijk + Mijk, n ni,knj. m , (1) 

where 

Lik = ]^ J qiU k rdv', (2) 



2 



4 ./,., 

and 



Nijk = t / qiUjU k r 2 dv', (3) 
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are elastic tensors. In the above expressions, 
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The derivatives of the interaction energy g(n, ft' , r) are evaluated on the reference state and the summation convention 
is assumed. Furthermore, u = r/r, and hence x k — u k r. 

The integrations are performed over the interaction volume. Since the micelles we are considering are disk-like, 
we consider in the calculations an interaction volume in the shape of a disk. In assuming an interaction volume of 
this form we suppose that the interaction energy g(n,n',r) is different from zero in the region limited by two similar 
disks, whose inner part coincides with the micellar volume, and the outer part is defined by the long range part 
of the intermolecular interaction. For simplicity, the two disks are supposed to be similar, having the same shape 
anisotropy, R = 2a/d = 2A/D, where a and d refer to the inner (micellar) volume, whereas A and D refer to the 
outer volume. The bilayer thickness d is supposed to remain essentially independent of the temperature, according to 
the experimental results 

We will consider an intermolecular interaction in the mixed form |l2lll3| 

C 

g(n, n', r) — — ^ [n ■ ft' — 3e(n ■ u)(n' ■ u)] 2 , (6) 

where we have introduced a mixing parameter e such that for e = the interaction is reduced to the Maier-Saupe 
law[I3 

C 

3ms = — 7-(n-n') 2 , (7) 



whereas for e — 1 it is reduced to the Nehring-Saupe lawfl4| 



SNS = -^[n-n' -3{u-u')(u-n)} 2 . (8) 

In the above expressions we have assumed perfect nematic order, i.e., we assume that the scalar order parameter 
5 = 1. 

The details of the calculations of the elastic constants as functions of the eccentricity of the micellar (molecular) 
volume can be found in ref. 01 ■ To obtain the elastic constants as a function of the micellar shape anisotropy one has 
to perform the integrations over an interaction volume of discotic shape, characterized by a shape anisotropy R. To 
perform the integrations in order to evaluate the elements of the elastic tensors we choose the z axis of the coordinate 
system along ft. In cylindrical coordinates x — pcosO, y — psinO, and z, we have 

x pcosO y psm.9 z 
m = — = , u 2 = - = , us = -, (9) 



and, hence N = n ■ u = z/r, where r = \J p 2 + z 2 . In the calculations we need to evaluate integrals of the type [1? 



/= / f{p,9,z)pdpd6dz, (10) 

J v' 

with /(p, z) a continuous function of the variables. It is convenient to put the integrations in the form 
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f{p,6,z)pdp~ dz f{p,6,z)pdp\ , (11) 

J-d/2 JO J 

which is simply the integral over the outer volume minus the integral over the inner (micellar) volume of the same 
functions of ref. [Tflj . now in a cylindrical geometry. The integrations are easily performed and yield the following 
ratio for the elastic constants 



#33 _ P31 

#11 q 31 : 



(12) 



where 



p 31 = 4(i?(32(l + i? 2 ) 2 (l + 2i? 2 ) - 48e(l + 8# 2 + 1LR 4 + 4R 6 ) 

+ 3e 2 (3 + 136i? 2 + 141i? 4 + 48i? 6 )) + (32 - 48e + 9e 2 )(l + R 2 ) 3 arctan(l/i?)) 

<73i = (i?(64(l + i? 2 ) 2 (3 + 2i? 2 ) - 128e(3 + 8i? 2 + 9i? 4 + 4i? 6 ) 

+ 9e 2 (21 + 56i? 2 + 123i? 4 + 48i? 6 )) + 3(64 - 128e + 63e 2 )(l + R 2 ) 3 arctan(l/i?)). (13) 



Let us consider first the case e = 0, which corresponds to the Maier-Saupe interaction. In this limit we obtain 

#33 _ 2(i? + 2R 3 + (1 + R 2 ) arctan(l/i?) 
#il ~ i?(3 + 2i? 2 ) + 3(1 + R 2 ) arctan(l/i?) ' 

In the other limit, for e = 1, which corresponds to the Nehring-Saupe interaction, we obtain 



(14) 



#33 _ 4(7i? - 152i? 3 - 55R 5 - 16R 7 + (7 + 21i? 2 + 21i? 4 + 7i? 6 ) arctan(l/i?) 

#Ti ~ 3i? + 8R 3 - 403i? 5 - 48i? 7 + (3 + 9i? 2 + 9i? 4 + 3i? 6 ) arctan(l/i?) ' ( ' 

These results indicate that the ratio of the elastic constants is an increasing function of the micellar shape anisotropyjfj 
but it tends to a saturation in both cases (see Fig.^). In fact, when we consider R — > oo we obtain 

#33 = 4(2 -3 £ ) 2 

#11 8-32e + 27e 2 ' V ' 

and the saturation value obviously depends on the mixing parameter e . For e = 0, #33/#n = 2, whereas for e = 1, 
## 33 /#n = 4/3. Of course very high values for R make no physical sense. However, as we will discuss below, in 
order to account for the experimental values of the ratio #33/#n in a discotic phase, the saturation value obtained 
in the above calculation cannot be restricted to the interval 4/3 < #33/#n < 2. To account for these experimental 
data we consider the results for the interaction law having a mixing parameter e = 0.35. For this value, the predicted 
saturation value for the ratio #33/#n rj 33.6. It is important to stress again that the "pure" Maier-Saupe interaction 
and the "pure" Nehring-Saupe interaction cannot furnish a ratio of this magnitude, since, as stressed before, they 
have saturation values which are small, when compared with the ones experimentally found. 

Our procedure to obtain the trend with temperature for the ratio #33/#n is the following. We assume that the 
dependence of this ratio with the micellar shape anisotropy is given by Eqs. (|12f> and (|13|) . From the experimental 
data of ref. § we obtain the micellar shape anisotropy R for the corresponding value of the reduced temperature 
T/Tni, where Tni is the nematic-isotropic transition temperature. This result in turn furnishes the trend of K^fKw 
versus the reduced temperature as predicted by the theoretical calculation we have performed. 

In Q the splay and bend elastic constants have been measured, as a function of the temperature, for the mixture 
DAC'l, NH4CI, and water, in a nematic discotic phase. From these experimental data it is possible to obtain the ratio 
#33/#n, studied by the authors, as a function of the temperature for two compositions: (mixture 1) DACl 7.00, 
NH4CI 2.53, H2O 90.40 mole %; this sample presents the following sequence of phases: neat soap, 35 °C, nematic, 55 
°C, isotropic; (mixture 2) DACl 7.57, NH 4 C1 2.73, H 2 89.70 mole %; neat soap, 45.5 °C, nematic, 65.5 °C, isotropic. 
Their results shown that, in the vicinity of the neat soap phase the bend elastic constant diverges while the splay 
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constant shows no critical behavior. Among the measured values for the ratio, the higher one is K33/K11 ~ 5.3 for 
mixture 1 and K33/K11 w 19.3 for mixture 2. 

The ratio K33/K11 versus the reduced temperature of the theoretical calculations (for a mixing parameter e = 0.35) 
and the experimental data for mixture 1 are shown in Fig. The ratios, both experimental and theoretical, fall on 
the same curve. A similar behavior has been observed also for mixture 2. With this procedure we have obtained, 
in a phcnomcnological way, the dependence of the micellar shape anisotropy R on the reduced temperature for a 
discotic nematic phase. Furthermore, X-ray diffraction results 6j performed in a nematic discotic lyotropic mixture 
(potassium laurate - KL, 1-decanol-DeOH and D 2 - mixture 3) show that the temperature variations of the shape 
anisotropy and of the orientational order parameter are similar within the experimental error. The bilayer thickness 
is found to remain independent of temperature and in this lyotropic system the nematic discotic exhibits a reentrant 
behavior. 

The connection between experimental and theoretical results for the ratios R versus reduced temperature are 
presented in Fig. El Note that, on coming from the discotic nematic to the isotropic phase our results agree well with 
earlier measurements of the shape anisotropy on the mixture 3, particularly in vicinity of ncmatic-isotropic phase 
transition. A more precise comparison between our theoretical calculations with the experimental data is difficult 
since there are a few data available on the measurements of the shape anisotropy in lyotropic systems. The results show 
another interesting physical behavior since the temperature dependence of the ratio comes from the micellar shape 
anisotropy dependence on the temperature instead of the scalar order parameter dependence. In fact, the temperature 
dependence of the bulk elastic constants comes (approximated) from its quadratic dependence on the scalar order 
parameter. Therefore, this dependence is not decisive in the ratio of two bulk elastic constants of thermotropic NLC 
(it is particularly weak near the nematic-isotropic transition 17]). But the experimental data reported above for 
lyotropic NLC shows that this dependence exists and is also very strong. In this manner, it seems to be meaningful 
to consider that the temperature dependence for lyotropic NLC can be due mainly to the micellar shape anisotropy 
behavior with the temperature. Another important point to be considered here is the introduction of the mixing 
parameter e on the intermolecular interaction law among micelles of the lyotropic medium which seems to be the key 
of this study. The choice of e = 0.35 is just an indication that "pure" laws like the Maier-Saupe and Nehring-Saupe 
ones can not account for the experimental results, and the mixed interaction has to be taken into account. 
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FIG. 2: The ratio of bend and splay elastic constant vs. temperature for mixture 1 as predicted by this work (x). The 
experimental data (A) comes from ref.||. Tni is the nematic-isotropic transition temperature 
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FIG. 3: The trend of the micellar shape anisotropy R vs. temperature. The curves for mixtures 1 (A) and 2 (□) comes from 
our theoretical predicitions (best fits) whereas the curve for mixture 3 (0) is the experimental results from ref.|g]. 
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